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STABILITY OF STEADY STATES OF THE NAVIER-STOKES-POISSON 
EQUATIONS WITH NON-FLAT DOPING PROFILE 

ZHONG TAN, YANJIN WANG, AND YONG WANG 


Abstract. We consider the stability of the steady state of the compressible Navier-Stokes- 
Poisson equations with the non-flat doping profile. We prove the global existence of classical 
solutions near the steady state for the large doping profile. For the small doping profile, we 
prove the time decay rates of the solution provided that the initial perturbation belongs to 
with 1 ^ p < 3/2. 


1. Introduction 

The dynamics of charged particles of one carrier type (e.g., electrons) can be described by 
the compressible Navier-Stokes-Poisson equations: 

'dtp + div(pu) = 0, 

dt{pu) + div(/9n (8) n) + Vp{p) — pAu — {p + p)V div u = pVc/), , , 

A(j) = p — b, 

Sp,u) |t=o= {po,uo). 

Here p = p{t, x),u = u(t, x) represent the density and velocity functions of the electrons respec¬ 
tively, at time t ^ 0 and position x G The pressure p = p{p) is a smooth function with 
p'{p) > 0 for p > 0. We assume that the constant viscosity coefficients p and p' satisfy the 
usual physical conditions 

p>0, p'+‘^p^0. (1.2) 

The self-consistent electric potential (j) = 4>{t, x) is coupled with the density through the Poisson 
equation, where the fnnction b = b[x) is the doping profile for the ions. We assume that 6 is a 
smooth function satisfying 

b{x) > 0, lim b{x) = 6 > 0. (1.3) 

|a:|—)-+oo 

For the pressure law p[p) = p^ with the adiabatic exponent 7 > 3/2, the global existence 
of weak solutions was obtained by [8] when the spatial dimension is three in the framework of 
Lions-Feireisl for the compressible Navier-Stokes equations HIE]. This result was later extended 
by dH to the case 7 > 1 when the dimension is two, where the anthors introduced an idea to 
overcome the new difhcnlty caused by that the Poisson term pV(p may not be integrable when 
7 is close to one. The large-time behavior of weak solntions towards the steady state was also 
considered in m- When the doping profile is flat, i.e., b{x) = 6, the steady state of (jl.ip is 
the trivial constant one (p, 0,0) with p = b. The uniqne global solntion aronnd this constant 
state in was proved by [lOj in the framework of Matsnmura-Nishida for the Navier-Stokes 
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equations m- Moreover, when the initial perturbation {po — p, uq) is small in L^, the decay 
of the solution to (11.11) was also obtained in m- 

II(p-p)WIIl 2 ^ (1 + ^)"^ and ||u(f)||^2 < (1 + t)“i (1.4) 

This may imply that the presence of the electric field slows down the time decay rate of the 
velocity with the factor 1/2 compared to the Navier-Stokes equations [T 2 l[T 3 ]. It was proved by 
transforming the system o into the Navier-Stokes equations with a non-local force 


dtp + diY{pu) = 0 

dt{pu) + diY{pu + Vp(p) — pAu — {p + p')^ divrt = pVA~^[p — p). 


However, the author in gave a different (contrary) comprehension of the effect of the electric 
field on the time decay rates of the solution. When the initial perturbation (po — V(/>o) 

belongs to with 1 < p ^ 2, the decay of the solution to (II.1|) was obtained in [IB] : 




2 and \\u{t)\\^2 < (1 + t) 



( 1 . 6 ) 


In this sense, the electric field enhances the time decay rate of the density with the factor 1 / 2 ! 
This can be understood well from the physical point of view since we get an additional dispersive 
effect from the repulsive electric force. 

In this paper, we will study the asymptotic stability of the steady state of the system 
with the non-flat doping profile b{x). A steady state {ps, (ps) with ri* = 0 of (II.ip must satisfy 


I Vp{ps) = Ps'^fps, ^ 

[Afps = p^-b. 

We will record the existence and uniqueness of the solution to (II. 7 p in Proposition 12.11 

Notations. We use TP(M^), I ^ p ^ oo to denote the spaces with norm and IT^’P(M^) 

to denote the usual Sobolev spaces with norm ||•||^yfc,p, and with £ S M stands 

for the usual spatial derivatives of order we allow that I <d oi i \s not a positive integer. 

Throughout this paper, we let C denote the universal positive constants. We will use A < B 
if A ^ CB and A > H if A ^ CB, and we may write ^A + B < D for ^A -|- CB < D. For 
simplicity, we write ||(A,H)||^ := ||A||^ -h ||H||^ and // := f^s fdx. 

Our first main result of the global solutions to dnn near the steady state for the large doping 
profile is stated as the following theorem. 

Theorem 1.1. Assume that Vb{x) G with k ^ 3, and {ps,(ps) of (11.71) is constructed in 
Provosition \2.1[ If ||(po — Ps,uo)\\fjk + ||V“^(po “ /^^f)||L 2 sufficiently small, then there exists 
a unique global solution {p,u,'V(p) to the system (II.ip such that for all t^O, 

||(p - Ps,u){t)\\]jk + ||(V(/> - Vcps){t)\\l 2 + (^IKp - Ps){T)\\]jk + llV'u(r)ll^fe) dr 

^ C {\\{P0 - Ps,Uo)\\Hk + ||V“^(po -p^)||^2) • (1-8) 

Remark 1.2. In Theorem \1.1[ the initial condition on the smallness of ||V“^(po — p<i)||j ;^2 is 
required so that the initial potential energy ||V(/>o — V(ps\\L'^ is small. Indeed, since (p—(ps satisfies 
the Poisson equation A{(p — (ps) = P — Ps, we have 

\\V(po - VcPsWl^ = ||VA"^(po - Ps)\\l2 = ||V"^(po - Ps)\\l2 ■ 

Such condition can be guaranteed by that, for instance, ||po “ PsH^e/s is small. 

Our second main result of the time decay rates of the solution to dni) towards the steady 
state for the small doping profile is stated as the next theorem. 
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Theorem 1.3. Assume that the assumptions of Theorem \1.1\ hold for k ^ A. If ||V 6 ||j|^fc + 
||6 — b\\^r with 1 < r < 3/2 is sufficiently small and ||(V“^(/3o ~ with 1 ^ p < 3/2 

is finite, then for 0 ^ ^ 1/2, 


V\p-ps)it) 


—l( 1 i\ ^ 1 

^ 2 Vmax{p,r} 2^ 2 2 


and for 0 ^ ^ 3/2, 


and 


V^u(i) 


Hk-t 


^ Co{l + 1) 


3 ( _ 1 _ 

2 y max{p,r} 


£ 

2 


\\{p-Ps,u){t)\\^^ i^Co{l + t) 2 


1 l\ 

max{p,r} 2 J 


3 

4 


(1.9) 

( 1 . 10 ) 

( 1 . 11 ) 


where Cq is a positive constant depending on the initial data. 

Remark 1.4. In Theorem ll.A the initial condition on the boundedness of \\'\7~^{po — ps)\\j^p 

with 1 ^ p < 3/2 is required. If p > 1, then such condition can be guaranteed by that, for 

instance, po — Ps = div/ for some f G L^, thanks to the singular integral theory [Si- 

Remark 1.5. Theorems 11.11 a,nd, \l . A extend, the previous results of US]. Note that the argument 
of proving the time decay of the solution in [16], which follows a pure energy method introduced 
in |5|, highly depends on that the doping profile is flat. To show the time decay of the solution 
for the non-flat doping profile, we shall need to employ a different argument as explained below. 

Remark 1.6. In mm, the authors proved the time decay of the compressible Navier-Stokes 

equations with an external potential force provided that the initial perturbation belongs to 

with 1 ^ p < 6/5. We may expect to employ our arguments of proving Theorem \1.A to extend 
the range of p therein to be 1 ^ p < 3/2. The key point is to introduce the fractional derivatives 
in the study of the time decay as already seen from (|1.9p - (|l.inp . 

Theorems o and oi will be proved in Sections [SHU respectively. To prove Theorem ll.il we 
will reformulate the system dm) into (|3.5p for the perturbation Q = p — Ps and u = u. To derive 
the energy estimates, the difficulty is caused by the terms on the left-hand side of (13.5|) . More 
precisely, we can not directly control as in [S] the terms resulting from when the differential 
operator V* commutates with the functions of ps since ps may not be close to a positive constant 
in the current case. To overcome this difficulty, we first notice that these commutator terms do 
not appear when 1 = 0, which allows us to derive the zero-order energy estimates as stated in 
Lemma 13.11 When I = 1,2,... ,k, we will carry out the delicate analysis so that we can control 
these commutator terms with an (large) error term as stated in Lemma 13.21 briefly speaking, 
I|v«IL 2 . However, this error has been controlled by the previous step. Hence, after recovering 
the dissipation estimates of g by Lemmas I3.3H3.41 we can close the energy estimates. 

To prove Theorem 11.31 we will reformulate (13.51) into the system (14.ip with constant coeffi¬ 
cients. Since ps is close to p, we can improve the energy estimates in the proof of Theorem 1 1.1 1 to 
deduce the estimates (|4.59l) . which implies that the decay of ||-|- ||V^u(f)||^fc_£ 

can be obtained from the decay of ||V^“^p(f )||^2 + ||^^'“( 0 ||l 2 + 11(^1)'^)(^)llioo for 0 ^ ^ 3/2. 

On the other hand, using the linear decay estimates of the linear Navier-Stokes-Poisson equa¬ 
tions with constant coefficients, we can derive the estimates conversely. This interplay would 
then be closed by the smallness of the solution and the doping profile. Finally, we may remark 
that if we do not introduce the fractional derivatives, then we can only choose £ = 0,1 in p4.59p . 
which would result that we could only prove Theorem 11.31 for 1 ^ p < 6/5. 

The rest of this paper is organized as follows. In Section [21 we will prove the existence of 
the stationary solution to dEZD. We will prove Theorem 11.11 and Theorem 11.31 in Section [3] and 
Section 01 respectively. Some analytic tools will be collected in Appendix. 
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2. Steady state 

In this section, we record the following existence and uniqueness of the solutions to ()1.7p . 

Proposition 2.1. Assume that b{x) is a smooth function satisfying ()1.3I) . Then there exists a 
unique classical solution {ps,4>s) to ()1.7p . Moreover, 

• ps has the positive upper and lower bounds, i.e., 


0 < inf b{x) ^ Psix) ^ sup 6(x) < oo; 




cGR3 


( 2 . 1 ) 


ifVb G with k ^ 3, then there exists a constant C depending on ||V6||j:^fe such that 
sufficiently small, then 

if further ||6 — < oo with 1 < r < oo, then 

\\ps-p\\w^.r < \\b-b\ 

where p = b. 


IL’’ ’ 


( 2 . 2 ) 

(2.3) 

(2.4) 


Proof. The existence and uniqueness of the classical solutions to (|1.7p satisfying the first three 
assertions were proved in [6] for k = 3, but the case k ^ 4 can be handled in the same way and 
so we omit the proof. We may then focus on proving the last assertion a priori. 

To this end, setting h'{s) = pfs)/s, we deduce from (11.71) that 


div ih'ips)Vps) = ps-b. 

Writing f = ps — p with p = b, we rewrite (12.5p as 

-h'{p)Af + / = div {{h'ips) - h'ip)) V/) +6-6. 

It then follows from the standard elliptic theory m on (ESP that for 1 < r < oo, 

\\f\\w^.r < ||div iih'ips) - h'ip)) V/) + 6 - 6||^. 
<\\^Ps\\H>‘\\f\\w^,r + \\b-b\\^^. 

This implies (12.4p since is small by ()2.3I) . 

3. Global solution with large doping profile 


(2.5) 


( 2 . 6 ) 


(2.7) 

□ 


In this section, we will construct the global solutions near the steady state to CIP for the 
large doping profile. For this, we define the perturbation by 

Q = p- Ps, u = u, ^ = (j)- 4 >s- (3.1) 

In order to reformulate the problem ESP properly, we introduce the enthalpy function 

p'is) 


h{z) = f 


ds. 


We also introduce the Taylor expansion 

He + Ps) = Hps) + h'ips)e + TZ, 
where the remainder TZ is given by 


(3.2) 

(3.3) 
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Then the problem (11.111 can be reformulated into the perturbed form of 

/ 

dtQ + divipsu) = - div{Qu), 

dtu + V{h'( ps)q) - (pAu + (p + p')y divu) — V<h 

Ps 

^ =—u-Vu— 'S/TZ+{ -) (uAu + (u + n'lVdivtt) , 

\Q + Ps PsJ 


A^ = q, 

{Q,u) |t= 0 = 

3.1. Energy estimates. In this subsection, we will derive the a priori estimates for the solu¬ 
tions to the Navier-Stokes-Poisson equations (j3.5p by assuming that for sufficiently small <5 > 0, 

||(^>,«)(t)||^. + ||V«h(t)||^. ^<1. (3.6) 

We first derive the zero-order energy estimates for the solution itself. 


Lemma 3.1. It holds that 
d 


dt 


{h'{p,)e^ + Ps\u\^ + |V<I>|2) + \\Vu\\l, < 6 WeWl^. 


(3.7) 


Proof. Multiplying the first two equations in (|3.5p by h'{ps)g and psU respectively, summing up 
them and then integrating over M^, we obtain 

1 d 


2 dt 


{h\ps)Q^ + ps\u\^) + p ||Vu||^2 + {p + p') ||divu||^2 - J PsU- V$ 

= - J {h'{ps)Qdiv{Qu) + PsU ■ {u ■ Vtt) -I- PsU ■ viz) 

+ I PsU ■ ( — - - ^ ) (pAu + {p + p')V div u) . 

J \Q + Ps PsJ ^ 

Here we have used the integration by parts to have the cancelation: 

/ {h'{ps)Qdiv(psu) + PsU ■ S/{h'{ps)Q)) = 0 . 


(3.8) 


(3.9) 


To estimate the Poisson term on the left-hand side of (13.8p . we integrate by parts by several 
times, use the first equation and the Poisson equation in (13.5p . by Holder’s and Sobolev’s 
inequalities and the a priori bound (13.6p . and employ the standard elliptic estimates ||V ^$||^2 ^ 
II for A: ^ 2, to deduce 

— J PsU ■ V$ = J div(/9si^)‘h = — J dtp^ — J div(^3tt)<l> 

= - J dtA^^ + j 6^-^^ = J + J gu-V^ 




1 d 
2dt 


iv^r-ci 


Il 2 \W\\l3 IIv^IIl 6 




li- 


We now estimate the terms on the right-hand side of 
inequalities, (13.611 and (|2.ip . we obtain 

- / {h'{ps)gdiv{gu) + PsU ■ {u-Vu)) 


(3.10) 

By Holder’s and Sobolev’s 


/ 


+ PsU- 


1 


1 


j [pAu + {p + p')V div u) 


.Q + Ps Ps J 
^ WpWl^ l|Vt-||i3 ||u||i6 + \\g\\L 2 WpWloo ||divu ||^2 
+ \W\\l<^ II^IIls II^“IIl2 -I-11^11^6 \\q\\l2 ||v^m||j^3 
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< 


5(iieiii2 + iivn|ii.). 

For the remaining term involved with V7^, we integrate by parts to have 

— J psU ■VTZ = J V ps ■ uTZ + ps div uTZ 


(3.11) 


^ II^IIloo ||£>||i 2 + ||divu ||^2 \\q\\l 2 ll^-llio 


< 6 


1,2 


+ l|Vu| 


1,2 


(3.12) 


Here we have used the fact from (13.3^ - 113.4p that TZ = 0{q'^). 

Plugging the estimates (I3.10p - (j3.12p into ()3.8p . since 6 is small, we then conclude (13.7p . □ 

We then derive the energy estimates for the derivatives of the solution. 

Lemma 3.2. For I = 1,... ,k, we have that for any e > 0, 


dt 


I (^h'{p,)\VQ\^+Ps\Vuf) + 




+ \\Q\fL2 ) + Cs 


L 2 

V'u 


L 2 


+ iivu|ii 2 


(3.13) 


Proof. Applying V* to the first two equations in (|3.5I) and then multiplying the resulting iden¬ 
tities by h'{ps)'S/^g and psV^u respectively, summing up them and then integrating over we 
obtain 

ijt / {h'iPs)N'Q\^+Ps\V^u\^)+ j [h'ips)VQV^diY{p,u)+PsVu-V^Vih'{ps)Q)) 

— J ps^^u ■ ^— (/rAti + {p + p)V div u)^ — J psV^u ■ V^V$ 

= - j h'{ps)V^QV^ div(^n) - j psV^u ■ V\u -Vu)- j psV^u ■ V'V7^ 


+ 


jp.V 


1 


u • ( (-) (pAu + (p + p')V div u) ) . (3-14) 

Wq + Ps Ps J J 

First, we estimate the terms on the left-hand side of (|3.14l) . By integrating by parts by 
several times and employing the commutator notation (IA.2p . we have 


J {h'{ps)V^QV^ d\v{psu) + psV^u • vV(h'(p, 

= J (^h'{ps)S/’-gV’' div{psu) - div{psS/’-u)V’‘{h'{ps)g)^ 

= J (h'(p,)V'^divQv^p,] u)V [v^/^'(p,)] 


= / h'{p,)Vg(\v\ps 


divu + 


+ / p.V'u- ( V^,h'{ps) 


Vg + 


v‘,Vp, 

V\Vh'ips) 


(3.15) 


For the second term on the right-hand side of (I3.15p . we employ the commutator estimates 
(IA.3I) to have: for / = 1, 


J h'{ps)Vg [V, Vps] u < \\Vg\\^, || [V, Vp,] u||^2 
< l|Vp|L2 lIVV.ILs hhe < ||Vp||i2 ||Vu||^2 


(3.16) 


for I ^ 2, 


h'{ps)V^g V^Vp, 


u 


< 




L 2 


v^Vp, 


L 2 
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< 


< 

/'N-/ 


v'e 

Vg 


L 2 


L 2 




L3 




1,6 


+ 


v'+Vs 


L 2 


ru 


Vu 


L2 


+ lIVnii 




< 

rsj 


V'^? 


L 2 


V'u 


L 2 


+ liVnIli. 


Note that we have used (j2.2p and the interpolation from Lemma lA.ll ||V^n ||^2 ^ 
||Vn ||£2 for I ^ 2. We may then conclude that this term can be bounded by 


Vg 


L 2 


Vu 


1,2 


+ l|Vu| 


L 2 


We can derive the same bound for the first term that for 1 = 1, 

J h'{ps)Vg[V,ps]dwu < \\Vg\\^2 ||[V,ps] divn||^2 
< l|V^||i 2 llVp.ll^^ ||Vn ||^2 < ||V 0||^2 ||Vn ||^2 ; 

for I ^ 2 , 

/ hf(o^)\/^o VKn. div?/, < VK n. div?y, 

L 2 


1 h'{ps)Vg 

^\ps 

divn < 

V'i? 

L2 

V^p. 

divn 


< 


< 


^2 (ll^Psllioo Mivn 




L 2 




1,2 


L 2 

< 


+ 


VVs ^Jldivull^a) 


V'i? 


1,2 


V'u 


L 2 


+ 


Similarly, for the third and fourth terms, we have: for 1 = 1, 

J psVu-{[V,h'{ps)] Vg+ [V,Vh'(p,)] g) 

< ||Vu||i 2 (||[V,h'(p,)] Vg\\^, + ||[V,Vh'(p,)] g\\^,) 


< 


for I > 2 , 


||Vu ||^2 (||Vh'(p,)||^^ \\Vg\y + \\V^h'{p,)\\^, ||£,||^e) < ||Vu ||^2 ||V ^>||^2 ; 


PsVu-(\v\h'{ps) 


Vg + 


< 

rs_/ 


< 

rs_/ 






V',h'(p.) 


V'i? 


+ 




< 


Vu 


1,2 


(||Vh'(p,)| 

,2 (I|v'^'(^«)IL3 


V',Vh'(p,) 

Vu 

+ 


L 2 

I u' 


v^v/^'(p,) 


L 2 


||V^y||^3 




< 


+ 

L 6 

V'u 


V^+^h'(p, 


1,2 


llelb 


L 2 




Hence, we may bound the third and fourth terms by 


Vu 

{ 


L2 V 


Vg 


L 2 


+ I|V£'IIl2 


We may thus conclude that 

J {h'{ps)VgV d\v{psu) + psVu ■ VV{h\p, 


> - 


+ ||V^?||^ 2 ) ( 


Vu 


L 2 


+ 


Next, we compute the following term 


u-V 


1 

Ps 


/ 


V divM ) = — / Ps^u • V‘ ( V ( — divM ) — V 


/.,v 

= j (^div (p,V'n) V' div+ p,V'u • V' div: 

= J |V^ divn|^ + J psV divn 


div u 


divn 


(3.17) 

U\\^2 + 


(3.18) 


(3.19) 


(3.20) 


(3.21) 


(3.22) 
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+ / Vps • V 


divu^ + J psV’'u-V’' 


(3.23) 


For the second term on the right-hand side of (j3.23jl . we employ the commutator estimates 
(1A.3P to obtain: for / = 1, 


PsV divu 

< llv^-" 


V,- 

Ps. 


divu < IIV divu ||^2 


[v,- 

divtt 

L ps\ 



L 2 


L 2 


V 


L°° 


l|divu||^2 < ||V\||^2 l|Vu||^2; 


for I ^ 2 , 


PsV divtt 


V^ — 

divu < 

div u 


v^ — 

div u 

Ps. 



1,2 

Ps. 



1,2 


< 

rs-/ 


< 







V 

V'u 


1 

Ps 

L 2 


+ llVull 


m 


Miv 

u 

+ 

v' 




L 2 

\ 

Kpsj 

) - 


( 

Vu 

J 



L 2 V 



L 6 


||divtt| 


L3 


L 2 


+ l|Vn ||^2 


For the third term, we employ the product estimates (IA.4I) to obtain: for I = 1 

J Vps • VuV ^— divu^ < l|Vu ||^2 V ^— divu 

||V divM ||^2 + 


< 


< 


l|Vu||^2 


V( - 

Ps 


L 2 

||divu ||^6 


|V^u||^2 ||Vu||^2; 


for I ^ 2 , 


Vu-V f — divu ) < 

Vu 


v' 

f— divtt 1 

\Ps J 


L 2 


\Ps J 


< 

r\-i 


< 


Vu 

( 


1,2 V 

Vu 



1,2 V 


1 


V div u 

+ 

V 

f-) 

Ps 

L°° 


L 2 


\PsJ 


L 6 


L 2 

||divu||^3 




L 2 


+ 


Similarly, for the fourth term, we have: for / = 1, 

I psVu • V (v (^1) divu) < ||Vu ||^2 


< 


< 

rsj 


llVull 


L 2 


V( - 

Ps 


|V divu||j ;^2 + 


V ( V ( — 1 divM 

Ps 
1 


L°° 


I — 

Ps 


1,3 


1,2 

||divu ||^6 


|v\||^2 l|Vu||^2; 


for I ^ 2 , 


j psV^u-V^ (^V (^^^divu^ < 




L 2 


VM V I — 1 divu 

Ps 


< 


< 


Vu 

( 


L 2 V 

Vu 

( 


L 2 V 




V div u 

+ 

V+^ f—) 

\PsJ 

L°° 


L 2 

\PsJ 


L 2 


L 2 

||divu||^o 




L 2 


+ llVull 


/f 2 


) < 

V^u 

( 




L 2 V 



L 2 


+ l|V^x||^2 


(3.24) 


(3.25) 


(3.26) 


(3.27) 


(3.28) 


(3.29) 


In light of (I3.24|) - (l3.29p and Cauchy’s inequality, we deduce that the last three terms in (I3.23P 
are bounded by that for any small constant p > 0 , 


P 




L 2 


+ Cn 


Vu 


L 2 
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We now compute the term 


1 


1 


Ps 


- I psVu ■ V‘ ( — Au ) = - / PsV'n • v' ( div ( —Vu ) - V ( — ) • Vu 


Ps 


Ps 


=y V (“) ■ 


= / 




V', 


Ps 


Vu 


+ / Vp, • V^M • v' ( — Vu ) + / p,v‘u • V‘ ( V ( — ) • Vu 


1 


Ps 


L. Til 


1 


Ps 


(3.30) 


So as for (j3.23p . the last three terms in (j3.3nD are also bounded by 


V'+^u 


L2 


+ Cn 


v'u 


L 2 


+ l|Vu| 


L 2 


Hence, we may conclude that 

— J psV^u • V* ^— (/rAu + (^ +/i')Vdivu)^ 


> 


V'+^u 

-c( 

v'u 


L2 V 



L 2 


+ l|Vu||i. 


We may simply bound the Poisson term as 


-j p,V'u-V'V$< 

Then, we obtain 


V'u 


L 2 


v'v$ 


L 2 


V'u 


- / p,V^u • V^V$ > - 


V'u 


L 2 


L2 


1,2 + 


L 2 


+ 






L2 


L2 


(3.31) 


(3.32) 


(3.33) 


Now, we estimate the nonlinear terms on the right-hand side of (|3.14p . By the commutator 
notation pA. 2 p . we have 


J h^ps)^^div(gu) = — J h'{ps)'V^QV^ {u ■ Vg + gdivu) 


j h'ips)^ 


u-VV^gV^g- / h'{ps)Vg 


V\u 


Vg- 


j h'{ps)V^gV^ {gdiv u). (3.34) 


By the integration by parts, we obtain 

- J h'{ps)u-VV^gV^g = -^ J h\p,)u ■ W\Vg\^ = ^ j div{h'{ps)u)\V^g\^ 


< 


|u||j;^oo -I- ||divu||^oo) 


V^g 


L2 


< (5 


By the commutator estimates (IA.3p . we obtain 


I h'{ps)V^ 


V\u 


< 




L 2 


< 5 


V^g 


1,2 




Vp< 

Vg 

2 

L2 


V'+^U 


V^g 

v'g 

L2 + 


L 2 


L 2 


(3.35) 


V\u 




1,2 


Vu 


L 6 


l|Vf?|| 


L3 


(3.36) 


By the product estimates (IA.4p . we obtain 

- J h'{ps)'^^gV^ {g div u) < 


Vg 


< 

rs_/ 




L 2 


L°° 


V^ div u 


L 2 


V (gdivu) 


Vg 


L 2 


L 2 

||divu||^oo 
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< 6 




L 2 


+ 




Hence, we deduce from (I3.34p - (I3.37I) that 

- j h'{ps)VQV^ div(H < <5 (ll V'i? 

We use the product estimates ()A.4p to obtain 


1,2 


+ 




L 2 


(3.37) 


(3.38) 


— / psV^u ■ ■ Vu) 


< 




< 


Vu 


L 2 


U\ 


< 5 


Vu 


L 2 


+ 


vVu 

2 

L 2 


L 2 


1,2 

+ 


V' (u • Vu) 


Vu 


L 2 


1,2 

l|Vu||^„ 




(3.39) 


By the integration by parts, we deduce 

— J ps^^u ■ VS/TZ = J Vps ■ VuVTZ + J psV^divnV^T^ 

II^PsIIls 


< 


V'u 


V'7^ 

+ 



v'Te 


L6 


L2 


L2 



L 2 


<5{ ||V£.||i 2 + 


L2 


+ 




L2 


(3.40) 


Here we have used the nonlinear estimates (jA.Sp of TZ stated in Lemma lA.31 
For the last term, we integrate by parts to have 


PsVu ■ V 


— -— 1 V^u 

Q Ps Ps 


L, . VT^-l 


= — VpsV^u ■ V 


1 


1 


Q + Ps Ps 


W - / PsV+V • V 


V7^—1 


1 


1 


Q + Ps Ps 


Vu 


< 

rsj 


Vu 


+ 


< 

rv-/ 


L 6 

7^+1,, 


V 


1,2 


1-1 


1 


1 


V' 


1-1 


Q + Ps Ps 

1 


-V^u 


L6/5 


1 


+ 

< <5 


V^V 
/+1 


L 2 


1 


Q~^ Ps Ps 

1 


Vu 


1,2 


Q + Ps Ps 

1 


L3 




1 


Q + Ps Ps 




L2 


+ 


V-^Q 


L2 


+ 


1,2 

V'+^u 
2 


+ 


1,2 


V 




1 


1 


L2 


+ 


V' 


1-1 


Q + Ps Ps 
1 


IV^ul 


L 2 


L3 


1 


S + Ps Ps 


\Vu\ 


1,6 


L3 


Hence, we deduce that 


/ 


Ps^u ■ V 
< (5 




1 


L 2 


(3.41) 


1 




+ Ps Ps 

2 

+ 


[pAu + (/i + /lO V div u) 


L2 




L2 


+ 




L2 


(3.42) 


Plugging the estimates (I3.22p . (I3.3ip - (|3.33p and (|3.38p - (|3.42l) into (13.141) . we then obtain 


/ h'(p,)iv'f?|2+p,iv^un + 


< 5 




L2 


+ 


V^u 


L 2 


+ 


V^+^u 
2 


2 

1,2 




L 2 


+ 


V^-^Q 


L 2 


+ l|Vp|li2 + llpilia 
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+ 




L 2 


l|Vn||2, + 




L 2 


+ 


VV 


L 2 


By the interpolation and Young’s inequality, since 6 is small, we deduce (I.S.l.Sp , 
Now we recover the dissipation for g, and we first deal with g itself. 


(3.43) 

□ 


Lemma 3.3. It holds that 

A 

dt 


u ■ V{h'{ps)g) + ||^>||i. + WVgf^, < + ||Vu||i.. 


(3.44) 


Proof. Multiplying the second equation in (13.51) by S/{h'{ps)g) and then integrating over by 
Cauchy’s inequality and the nonlinear estimates (jA.Sp . we obtain 

I dtu ■ V{h'{ps)g) + \\V{h'{ps)g)\\% - y' V4> • V{h'{ps)g) 

^ f||v2u||^, + ||u-Vu||^2 + ||V7e||^2+ ) ||V(/i'(p,)^>)|| 

V \Q + Ps PsJ ipj 


<\\vml,+6\\Vg\\l,. 


L 2 

(3.45) 


By the hrst equation in (13.5p , we integrate by parts for both t- and x-variables to have 
- J dtu ■V{h'{ps)g) =-^ J u-V{h'{ps)g) + j u ■V{h'{ps)dtg) 

u ■ ^/{h'{ps)g) + J h'(ps) div u{u ■ Vg + gdiv u + ps divu + u-Vps) 


fL 

dt 


'U\\l2 


^-Jtj u-V{h'{p,)g) + C\\Vi 
By the integration by parts and the Poisson equation in (13.51) . we obtain 

-j V^-V{h'ips)g) = I A^h'{ps)g = I h'{ps)g\ 

Plugging (I3.46|) - (|3.47p into (j3.45l) . since 5 is small, we deduce (j3.44l) by noticing that 

+ \\Vg\\l2 < J h'{ps)g^ + \\V{h'{ps)g)\\l ,. 

We thus conclude the lemma. 

We next derive the dissipation estimates for the derivatives of g. 


(3.46) 

(3.47) 

(3.48) 

□ 


Lemma 3.4. For I = 1,... ,k — 1, we have 


4 - f V'u-VV'£» + 

dt J 


S/^g 


1,2 


+ 




L 2 




L 2 


+ l|Vu||i. + ||^i||i.. 


(3.49) 


Proof. Applying to the second equation in (13.5p . multiplying the resulting identity by 
and then integrating over we obtain 


J Vdtu ■ VVg + 


I h'{p,)\v^+^ef 


- / v'v$ • 


< 

rs_/ 


L 2 


+ 




L2 


+ 


V {u • Vu) 


L 2 


+ 

V'+^7^ 

+ 


) 




L2 

\\g + ps PsJ ) 

1,2/ 



L 2 


(3.50) 


As (j3.46D . we deduce 

- J Vdtu • vVg = -^J + y V'u • S/S/^dtg 


dt 


V^u-W^g+ / divuV^ (u • V^3 + ^divu +/9<j divtt + u • V/9<j) 
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^ —— / V*U • VV*^ + 
at I 




L 2 


+ 


V\u ■ Vg) 


L 2 


+ 


V'(£»divti) 


L 2 


+ 


V^{ps divix^ 


L 2 


+ 


V\u ■ Vps) 


As (I3.47jl . we have 


- J ■VVq = j = 


Vg 


L 2 


L 2 


(3.51) 


(3.52) 


Plugging (|3.51l) - (l3.52p into (I3.50jl . and applying the product estimates ()A.4p . Lemma [A.II and 
the nonlinear estimates (ira as in Lemma Ea we obtain 


^ / V^M-VV'p+ / /i'(p,)|V'+^e|2 + 


2 

L 2 


< 

r\j 


+ 


V^+\h'{ps) 

L 2 


L 2 ■ 

+ II^^IIl2 + 


L 2 


By the commutator estimates (IA.3p . we obtain 


v'+Sh'(p,) 


< 


L 2 


|Vh'(p,)|| 


L 2 + 




V 


«+l 


L 2 


+ 


L 2 


V'+i/i'(p,) 


V 


«+l 


L 2 


(3.53) 


L3 


Il 6 


< 




L 2 


+ II^£'IIl 2 . 


(3.54) 


Plugging (I3.54P into (I3.53h . by Cauchy’s inequality and since 6 is small, we have 


4 - [ V^u-W^g + 

at J 




1,2 


+ 




L 2 


< 




L 2 


+ ll^^llia + 




L 2 


+ II^£'IIl 2 • 


By the interpolation and Young’s inequality, we further deduce (I3.49P from (13.551) . 


(3.55) 

□ 


3.2. Proof of Theorem II.IL In this subsection, we will prove Theorem ll.il Multiplying the 
estimates p3.49p of Lemma 13.41 with / = A: — 1 by a small fixed constant ei, and then adding it 
to the estimates (I3.44p of Lemma [331 by the interpolation and Young’s inequality, we deduce 

2 


^ (/« • '^{h'{ps)Q) 


+ ei / S/’^-^u-V^g]+MU 


< 


V 


fc+i 


u 


L 2 


+ llVulli^ . (3.56) 


Multiplying (13.561) by a small fixed constant 62 and then adding it to the estimates p3.13p of 
Lemma 13.21 with I = k, we obtain 


^(^1 h'{ps)\\/^g\‘^ + Ps\V’^u\^ + e 2 (I u-V{h'{p,)g) + 


ei / 


+ 


V 


k+l 


U 


L 2 


+ WqWhi^ ^ (^ + ^) IIs'IIh'' + 




L 2 


+ llVulli. . 


(3.57) 


Taking e sufficiently small in (I3.57p and since <5 is small, by the interpolation and Young’s 
inequality, we then have 


^ (^j h'{ps)\W^g\^ + + £2 (I n • V{h'{ps)g) + 


ei / V^-^u ■ V’^g 


+ l|Vn||l,. + ||f?||l,. < llVulli^. (3.58) 

Multiplying p3.58p by a small fixed constant 63 and then adding it to the estimates (13.7p of 
Lemma l3.11 since 5 is small, we deduce 

^ (/ + (^j h'{ps)\V^g\‘^ + Psl'^^ul'^ 

+€ 2(^1 u-Vih'{ps)g) + ei +\\Vu\\l, + \\g\\l,^0. (3.59) 
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Note that the expression under the time differentiation in (13.591) with properly small ei, €2 and 
€3 is equivalent to ||£»||^fc + + ||V‘h||^ 2 . Hence, integrating (13.591) directly in time, we 

obtain (HD. By a standard continuity argument, we then close the a priori estimates (13.61) if 
we assume at initial time that ||^)o 2^2 is sufficiently small. The global 

solution then follows by a standard continuity argument combined with the local existence of 
solutions. The proof of Theorem fTTl is thus completed. □ 


4. Time decay with small doping profile 


In this section, we will derive the time decay rates of the solution to HD towards the steady 
state. For this, we need to require that the doping profile is of small variation, that is, b{x) is 
near the constant b. Then according to Proposition 12.11 Psix) is near p with p = b. Owing to 
this fact, we may rewrite (133]) as: 


/ 

5t£i + pdivrt = -div((£i + /95 -p)u), 

dtu + h'{p)VQ — - (pAu + (p + p')S/ divtt) 

p 

= ^—-j---J i^pAu + {p + p')V div u) 

A4> = g, 

(q^u) |t=o= 


- V4> 


-u-Vu-vn- V{{h'{ps) - h'{p))Q), 


(4.1) 


4.1. Energy estimates. In this subsection, under the assumptions of Theorem 11.31 we will 
derive the further energy estimates for the global solutions to the Navier-Stokes-Poisson equa¬ 
tions (HD obtained in Theorem 11.11 By Theorem 11.11 the assumptions of Theorem 11.31 and 
Proposition 12.11 we have that 

II(£'!'“)(^)IIh'' + II^‘^(^)IIl 2 + IIps - /o|Ihj=+i ^ <5 (4.2) 

for some small constant 5 > 0 . 

We first derive the energy estimates for the derivatives of the solution of order I with I G [0, k]. 
Note that now I is not required to be an integer. 


Lemma 4.1. For 0 ^ I ^ k, we have 

^ I (h'(p)|V^£.|2 + /i|V'u|2 + |v'V<h|2) + 




L2 


< (5 




L2 


+ 


V^V4> 


L 2 


+ 


Vu 


+ llt'llL°° + 11^111°° + ll^^llia + ||V^t(||^3 j . (4.3) 


L 2 


Proof. Applying V* to the hrst two equations in ()4.1I) and then multiplying the resulting iden¬ 
tities by h'{p)'\/^g and pV^u respectively, summing up them and then integrating over M^, we 
obtain 

~ + pN^u\^) + j p\V+\\^ + {p + p')\Vd[vu\‘^ - j pV'u-VV4> 

= - j h'{p)VQV^ div((e + ps- p)u) - j pV^u ■ V\u • Vu) 

- j pV^u ■ VV7^ - j pV^u ■ V^V{{h'{ps) - h'{p))Q) 


+ / pV^u-V 


-- ) ipAu + ip + p')V div u) ) 

Q + Ps pj ^ V 


(4.4) 


First, we estimate the terms on the right-hand side of (|4.4I) . By the commutator notation 
(IA.2I) . we obtain 


J h'{p)V'‘QV^ di\{{Q + ps - p)u] 
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= - j h'{p)V'-QV’- {u ■ V{q + ps - p) + {g + ps - p) divtt) 


= - / h'{p)u-VVQVQ- / h'{p)VQ 




Vg 


- j h'{p)V^gV^ {u-V{ps-p))-j h'{p)V^gV^{{g + ps-p)d\Yu). 
By the integration by parts, we obtain 

- J h'{p)u-VVgVg = j h'{p)u ■V\Vg\‘^ = ^ j h'(p) dlv u\V^gl"^ 


< 

r\j 


lldivixll^o 
m 


V^g 


L 2 


< 5 


Vg 


L 2 


By the commutator estimates (IA.,3p . we obtain 

- [ h'{p)V^g\v^,u] -Vg^ Vg 


L 2 


S/\u 


< 


V'i? 


1,2 


L°° 


V^g 


L 2 


+ 




L 6 


Vp 

llv^lLs 


L 2 


< (5 




L 2 


+ 




L 2 


By the product estimates (IA.4D and ()4.2p . we obtain 

- j h\p)VgV (u • V(p, - p)) < II II {u ■ V(p, - p)) 


L2 


< 

/•N-/ 




L2 


11 ,°° 


< (5 


( 

V'p 

2 

+ 


V 


L2 



2 


L2 


+ 




L6 


I|v(p.-p)IL3 


1,2 


+ « 


and 

-f 


- I h'{p)V^gS/''{{g +ps - p)divu) < V'p V'((p + - p) divu) 


1,2 


< 

/•N-/ 




L2 




< (5 




1,2 


+ 




1,2 


div u 


+ l|V'“llz,3 


L 2 


+ 


V'p ||Vu||^oo + V\ps- p) 


L6 


Hence, we may conclude that 
- j h'{p)V^gV^ div((p + p* - p)u) < 6 


( 

V^p 

2 

+ 


V 


L2 



L 2 


+ lkllz,°° + II^'^IIls 


(4.5) 


(4.6) 


(4.7) 


(4.8) 


(4.9) 


(4.10) 


Next, by Holder’s and Sobolev’s inequalities, the product estimates (IA.4p and the interpola¬ 
tion estimates pA.ip . we obtain that for I = 0 , 

- I pV^u ■V^{u-Vu) = - j pu-{u- Vu) < ||u ||^3 llnll^a \\Vu\y < 5 || ; (4.11) 

V' {u ■ Vu) 


for / > 1 


- J pV^u • (tt • Vu) = 


< 


CO 

+1 

+1 

< 

I'N-/ 


LC 


V^u 


L 6 


L6/5 


L 2 






L 2 


L 2 


+ ll^i| 


I i+1 

11,2 




I 

i+1 

L 2 


, i+1 
V 2i u 


I 


1+1 

L 2 




1 

z +1 

L 2 


< (5 




1,2 


(4.12) 
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Hence, we may conclude that for I ^ 0, 

- [ pV^u • (u • Vu) < 5 




L 2 


(4.13) 


Integrating by parts and by the product estimates (IA.4p . the nonlinear estimates ()A. 6 p and 
(iia), we obtain 


- J pVu ■ Vvn = j pV'divnV'Te < 




V'77 


L 2 



1,2 




1,2 


+ 




L 2 


(4.14) 


and 


-1 pv^u-v^v{{h'{ps)-h'{p))e) = J pV^+^u-v\{h'ip,)-h'ip))Q) 

\s/\{h\p,)-h'{p))e) 

(\\h'{ps) - h'(p)\\^ 


< 

r\-/ 


< 




V^+^n 


L 2 


V\h'{ps)-h'{p)) 


< 6 


( 


2 

-|- 


V 


L2 



L 2 


+ 


(4.15) 


We now estimate the last term on the right-hand side of ()4.4p . For / = 0, we easily obtain 

[pu ■ ( -- ) (pAu + (u + u')V divu) 

J \Q + Ps pj ^ 

||V'u||^3<<^(l|Vn||i2 + ||v2u||^3). (4.16) 

O -r Oa O T 7 \ / 


< 

rs_> 


For / ^ 1, by the integration by parts, the product estimates (IA.4p and (14.21) . we obtain 

-- ) (pAu + (p + p')\/ div u) ) 

+ Ps pj V 

— [pAu + {p + p')V div 




Q + 

= - [ pV+\ ■ 

1-1 


Q + Ps 


< 

rs_/ 






L 2 


1,2 


V‘ 


1 


g + 

1 


-- ) (pAu + (p + p')V divtt) ) 

^Ps pJ ^ V 


Q + Ps P 




L 2 




V 


l-l 


1 


g + Ps p 


Iv^ul 


1,6 


L3 


< 5 




L 2 


+ V^u 


L3 


(4.17) 


Hence, we may conclude that for I ^ 0, 
'pV 


J pV^u ■ —-j- - [pAu + {p + p')V divtt)^ 


< 5 




L 2 " 


(4.18) 


Now, we turn to estimate the terms on the left-hand side of (14.41) . For the second term, we 
deduce from (na) that 


j /r|V^+^up + (/i + ^')|V^ div^ 




L 2 


(4.19) 


For the remaining Poisson term, we integrate by parts and use the first equation and the Poisson 
equation in ()4.ip to obtain 
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- j pV^u-V'V^> = j div(pix) 

= - j div((^ + ps- p)u) 

= - j - V\{q + Ps- p)u) • V'V^> 

|V'V^>|2+ [ V\{q + Ps-p)u)-V^V^. 


- \!L 

2 dt 

By the product estimates (IA.4P and ()4.2p . we obtain 

[ V\ips - p)u) ■ V^V4. < V\{ps - p)u) 


(4.20) 


L 2 


VV4> 


< 


(lbs -pIIj 


vb 


< 5 


VV4> 


L 2 


+ 


L 2 

Vu 


+ 


^\PS-P) 


L2 


\U\ 


L 2 

V'V4> 


L 2 


L 2 


+ b 


(4.21) 


and 


I vbH-vV$< ||vbH 

vb 


L 2 


V^V4> 


L 2 


< 


< (5 


Vb 


+ 


L 2 


V'V4> 


1,2 


V^V4> 


L 2 


+ 


Vb 


L 2 


+ 




L 2 


(4.22) 


( 

V'V4> 

2 

+ 

vb 

2 

+ 

vb 

V 


L2 


L2 



L 2 


+ b 


Hence, we may conclude that 

- J pV^u ■ VV4> I W4>p - C6 

(4.23) 

Consequently, plugging the estimates (I4.10p . (|4.13|) - (|4.15l) . (I4.18p - (l4.19p and (I4.23P into (|4.4I) . 
since 6 is small, we deduce (14.3|) . □ 

We now recover the dissipation estimates for g. 

Lemma 4.2. For 0 ^ I ^ k — 1, we have 


A. 

dt 


Vb-vvb + 


< 


V'+b 


L 2 


+ 


Vb 

V'+b 


L 2 


+ 


v'+b 


L 2 


L 2 


+ d (ibllioo + Iblli°° + ll^^llis + II. (4.24) 

Proof. Applying to the second equation in (14.11) . multiplying the resulting identity by g 
and then integrating over we obtain 


j vbiU-vV£» + j /i'(p)|v^+bP-y v'v$-v'v^ 


< 


V'+i {{h'ips) - h'{p))g) 


+ 




+ 


V' 


L 2 

1 


+ 


V'+b 


L 2 


+ 


v' {u • Vn) 


L 2 


- - Vb 


L 2 


.Q + Ps P 

As (j3.46p . we deduce 

- J Vdtu ■ vvb = -^ j 


v'+b 


L 2 


(4.25) 


dt 


vb • VVb + / VbivixVbii ■ V ^3 + ^divu + ps divu + u • Vp<j) 
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^ —— / V*U • VV*^ + 
at I 




L 2 


+ 


V\u ■ Vg) 


L 2 


+ 


V*(^3divtt) + V^(/9sdivix) + V\u-Vps) 


L 2 


L 2 


L 2 


(4.26) 


As (j3.5ip - ()3.53p . applying the product estimates ()A.4p and the nonlinear estimates ()A. 6 p as in 
Lemma l4.ll we obtain 


^ J v^«-vv'p + J /i'(p,)|V'+VP + 




1,2 


< 


( 


+ 


V 


L 2 



L2 




+(5 ||Vix||^3 + (5 II V^u ||^3 + 6 




L 2 




L 2 


(4.27) 

□ 


By Cauchy’s inequality and since S is small, we deduce (I4.24p from (j4.27p 

We now combine Lemmas I4.1H4.2I to derive the following proposition. 

Proposition 4.3. Let k ^ 3 and 0 ^ i ^ 3/2. Then there exists an energy functional 
equivalently to ||V^(^>,u, V4>)||^j,_£ such that 

d 


dt 


£^ + 


V^g 

2 

+ 




2 

+ 



Hk-e 


Hk-e \ 


L2 



2 

L 2 


\ l °° a II'^IIl° 


(4.28) 


Proof. Summing up the estimates (1331) of Lemma KT\ for from / = £ to A:, by the Poisson 
equation in (|4.1h . we obtain 


dt 


h'{p) V^g 


Hk-t 


+ P 


Vu 


Hk-e 


+ 


VW4> 


2 

Hk-e 


+ 




Hk-e 


< 6 


( 

V^g 

2 

+ 

VV4> 

2 




Hk-e 


L2 





(4.29) 


Summing up the estimates (14.2411 of Lemma 14.21 for from I = i to k — 1, we obtain 


A 

dt 


j V^u-VV^g + 

■1 1 'J 




vv 


Hk-e 


< 




Hk 


^ (ll£’lli°o + \\u\\Y + II^^IIls + ||^^'“IL3) • (4.30) 

Multiplying (I4.30p by a small constant e > 0 and then adding the resulting inequality to (14.291) , 
since S is small, we deduce that for 0 ^ ^ A: — 1, 


s I 


V^g 


Hk-e 


+ P 




Hk-e 


+ 


VV4> 


Hk-e 


+ e 




i<l€k-l 


+ 

Vg 

2 

+ 




Hk-e 



< 5 


VW^- 


L 2 


+ 




Hk-e 
2 


+ II£’IIl°° + + ll^^llia + ||V^u||^3 


L 2 


(4.31) 


We define to be the expression under the time derivative in (I4.3ip . Since e is small, £^ is 
equivalent to ||V^(g, u, V<f>)||^j,_£. Then we deduce that for 0 ^ ^ A: — 1, 


dt 


+ 


2 

+ 




Hk-e 



2 

Hk-e 


< 5 


V^V4> 


1,2 


+ 


V^u 


+ \\ q \\ l °° + ll^lli°° + ll^^llia + II^^^ILs 


1,2 


(4.32) 
































































































































































18 


ZHONG TAN, YANJIN WANG, AND YONG WANG 


Now we take 0 ^ i ^ 3/2. Then by the interpolation, we have 


llViill 


^ < 
L3 ~ 


V3/2 


U 


L 2 




1,2 


+ 




2 

L 2 


and 


;s 




< 


L 2 




L 2 


+ 




L 2 


Since 6 is small, (14.321) implies (I4.28p by (I4.33p - (l4.34p . 


(4.33) 

(4.34) 

□ 


4.2. Duhamel form analysis. In order to use the linear decay estimates for the linear system 
with constant coefficients, we will rewrite the Navier-Stokes-Poisson system dHD as the Navier- 
Stokes equations with a non-local self-consistent force in the following form: 

dtQ + pdivu = N^, 

dtu + h'{p)VQ — - [pAu + {p + p')y div u) — VA~^q = (4.35) 

, {q^ u) |i=o= (po,uo), 
where the “nonlinear” terms are given by 

= - div((p + ps - p)u) (4.36) 

and 


= —u ■ Vu — VTZ — Viih'ips) — h'{p))Q) -(- (-- ) ipAu -|- (^ -|- p')V divu) 

\q + Ps pj '■ 


(4.37) 

By the Duhamel principle, the solution [g, u) to the problem ()4.35l) can be expressed as 

{Q,u){t) = e~^^{Qo,uo) + [ (4.38) 

Jo 

Here the matrix differential operator A is defined by 


A = 


0 /odiv 

h'(p)V-VA"i -UpA + {p + p')Vdiv) 


(4.39) 


In light of the analysis in [TO] and [T 6 |, we have the followings about the time decay rates of 
the solution semigroup of the linearized system of (14.351) . 


Lemma 4.4. Let {g,u) = e ^^{go,uo). Then for l^p^2, q^2 and £ 0, we have 


and 


V^g 




Lt 


—3/ 

2\p q) 2 2 


(||(V ^go,uo)\\^p+ V\gQ,uo) 


Li 


Li 


_3/i_l __ 

<(l+t) Hp 1) 2 


(||(V ^£» 0 ,^io)||ip + V^(£>o,^io) 


Li 


(4.40) 

(4.41) 


Applying Lemma 14.41 to (14.381) . we obtain the following proposition. 
Proposition 4.5. It holds that for 1 ^ p,r ^ 2, q ^ 2 and £ ^ 0, 


_3 / 1_1 \_ 1_1 

V^g{t) <{l + t) Hp J 2 2 
t 


L9 


+ j {l + t-r) i) ^ ^ (^||(V-^iV\iV 2 )(r)||^,+ ||v^(iVi,iV 2 )(r)||^J dr 

° (4.42) 


and 


vVt)||^^ <(l + t) ^ (||(v-'^>o,^xo)|Lp + ||v^(^>o,^o) 


Li 
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+ 


[\i + t-Ty 

Jo 


3/l_l __ 

2\r qj 2 


(||(V-iiV\iV2)(r)||^, + ||v'(iV\iV2)(r) 


L‘1 


dr. 

(4.43) 


Finally, we record the following estimates of nonlinear terms appeared in (I4.42l) ~ (j4.43p . 
Lemma 4.6. It holds that for 1 < r ^ 2, 

\\{V-^N\Ny 


^ (<5 + IIPs - P\\w^,r) (I^IIl” + II^£'IIh 2 + ll^^'^lliLl) 


+ II£'IIl 2 


Q CJ 

V'* ru 


L 2 


+ II^IIl 2 


A ^ 

S/^-rU 


1,2 


forO^i^ 3/2, 


and 


V\N\N^) 1^^ < 5 (lliill^^ + \\Ve\\H2 + \\v\\\^,) ; 


||(^\ ^')ILoo < 5 (IIuIIloo + \\S/g\\H, + II 


(4.44) 

(4.45) 

(4.46) 


Proof. We will estimate the nonlinear terms term by term. First, for 1 < r ^ 2, by the singular 
integral theory [H], the identity (IA.9p in the proof of Lemma I A. 3 1 and Holder’s and Sobolev’s 
inequalities, we obtain 

||vMiv((£> + p^ - p)u)\\^, < ||(£» + Ps - p)u\\Lr 

" +\\ps-p\\Lr\\u\\L-°: (4-47) 


< ||^>||l2 ||u|| _i + IIp^ - p\\^r \\u\\loo < \\g\\L2 

L ^/r-l/2 


Q 6 

V'^ rU 


1,2 


I'w-Vu||^. < ||u||i 2 ||Vu|| < 


\u\\l^ 


V^--u 


L 2 


l|V7^||^. < \\h^'ips)gVg\\^^ + ||7^(/^') (V^. + Vp,)||^, 


< 11^11^,7^ IIVi?IL2 + II^'IIl- (W^qWl^ + W^psWl^) < s \\Vg\\^, ■ 


(4.48) 

(4.49) 


||V((h'(p,) - hyp))g)\\^^ < \\h\ps) - h'{p)\\^^ llVell^^ + ||V(/i'(p,) - h'(p))||^. ||^?||^^ 

~ IIt^s “ pIIvul^ ll^^/lliL 2 ; (4.50) 

1 1 ' 


e + Ps p 




< 

1 

1 

rs_> 

L'- 

Q + Ps 

p 


Iv^ul 


< 


L2 


'lV’— 1/2 


i||V' 


u\ 


L /1 • 


(4.51) 


These estimates (I4.47p - (j4.5ip give (j4.44p . 

Next, for 0 ^ ^ 3/2, by the product estimates (IA.4p . the nonlinear estimates (IA.5P and 

Sobolev’s inequality, we obtain 


A\Y{{g + Ps - p)u) <\\g + Ps- p\\l3 


< 


^ {\Mloo + ||v^u||^ 2 ); 


vyu • Vu) 

L 2 


< 


L2 ~ 



< (5 


L2 


1I^'“IIl6 + \\u\\]^3 



+ 


1,6 


< 


L 6 ~ 


V^^\g + Ps-P) 


L2 


IL°° 


^||V^ 


|JL 2 ’ 




L2 


< 


<5||Ve|lH2; 


S/^Vi{h'{ps)-h\p))g) 


1,2 


< 

/•N-/ 


V^^Hh'iPs) - h'ip)) Uhe + \\h\ps) - h'ip)l 


V" 


< 


1 


e + Ps p 


--W^u 


L2 




< 


1,2 


(4.52) 

(4.53) 

(4.54) 

6\\Vg\\^,-, (4.55) 




1 


II^^^IL 2 + 

P/ L°° 


1 

1 



Q + Ps 

p 

L°° 



.Q + Ps 

These estimates (I4.52p - (j4.56p yield (|4.45p . 

Note that the estimate (I4.46P can be obtained in a similar way. 


L2 


<s\\v^ 


H2 ■ 


(4.56) 


□ 
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4.3. Proof of Theorem II.3L In this subsection, we will prove Theorem 11.31 Let k ^ 3 and 
0 ^ f ^ 3/2. Adding ||V^(u, V<I >)||^2 to both sides of the estimates (j4.28p of Proposition 14.31 
we obtain 

2 

1,2 

2 


—yi! j. 
dt‘ + 


Vq 

2 

T 


2 

+ 

VV4> 


Hk-t 


jjk-e+i 



< 

rs_/ 


V"V4> 


L 2 


+ 




L 2 


II £*11 L°° ll^ll L°° 


Note that is equivalent to ||V^ ^£i||l 2 - Hence, we deduce from (j4.57ll that 






L 2 


+ 




L 2 


+ II£’IIl°° + II'“IIl° 


(4.57) 


(4.58) 


for some constant A > 0. By the Gronwall inequality, we obtain 




jjk+i-e 

JO 


+ 


V^u(t) 


Hk-l 


< 




Hk-e 


+ ||V$, 


oIIl 2 


-A(t-r) 


V'^-V(r) 


1,2 


+ 


V^u(r) 


L 2 


+ £• T 


+ \\u{t) 


dr. (4.59) 


We now prove (ll.9h - (ll.lip . So we let A: ^ 4, 1 < r < 3/2 and 1 ^ p < 3/2. For simplicity of 
notations, we denote 

^0 '■= ||(V“^po,^^o)||2,p + II(£'o,^^o)|Ih'= + II^^o|Il2 ) 

and we define 

1 1 ' 


C := ^ ( - 7 - 7 - ^ ) and 5 := h + \\ps - pWw^.r ■ 

2 \max{p, r| 2 / 

It turns out that we have to distinguish our arguments by the value of r. 
Case 1: 6/5^r<3/2. In this case, we define 


(4.60) 

(4.61) 




L 2 


+ 


+ ||0(f)||ioc + ||M(t)||LOO , 


M{t) := 


Hk-l/2 


+ 




Hk-3/2 ’ 


(4.62) 

(4.63) 


and 


J\f{t):= sup ((1+ r)^+4 (/:(r) + AI(r)) + (1+ r )^+2 ||^(r )||^2 + (1+t)^ ||u(r)||^ 2 ) . (4.64) 

O^TSjt A / 


(4.65) 


We take I = 3/2 in (|4.59l) to have, in view of (I4.62p - (l4.63p . 

M^{t) <e-^^Kl+ [ e-^(*-^)£2(r)dT. 

Jo 

We now estimate the time decay rates of C{t) by applying the linear decay estimates. By the 
estimates (I4.42P with ^ = 1/2 and g = 2 of Proposition 14.51 and using the nonlinear estimates 
(j4.44l) - (l4.45l) . in view of (I4.62I) - (I4.63I) . we obtain 

^ (1 + A) + [ {1 + t — + M-){t) dr 

Jo 


+ /(i + f-r)-i(^^H 

Jo 


Il 2 


V 


3 - 


^^ + \Ht)\\l 2 ru{T) dr. (4.66) 


By the interpolation, in view of p4.64p . we estimate 

73 — 


lL 2 


"u(r) +||u(r)||^2 "u(r) 


L 2 


< 

rv-/ 


2 -^ 


^-1 


V’^e(T) V^/2^(r) ||u(r) 


I --1 

I £2 


2 -^ 


/ + ||u(t )||^2 ru{T) 


L 2 


< 5 ^ ^(1 +t)( ^ 4)(r- ^(1 +t)( ^ 4 )^ r)j\ 7 (t)^ '•+(5(l+r) ^ 
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<(5(l + r) ^ (4.67) 

Here we have used the facts that d = (r — 2)/(4r — 4) ^ —1 and 4 — 3/r ^ 3/2 since r ^ 6/5. 
Hence, plugging the estimates (|4.67p into (14.661) . in view of (14.641) . we have 

ft 


^(1 + ^) + ^ + ^ — 4(l + r) ^ ^M{t)dT 

Jo 


< + (7^0 + 


(4.68) 


Here we have used the fact | (^ “ ^) + | > 1 since r < |. Similarly, by the estimates (I4.43p 
with £ = 3/2 and q = 2, and (14.42l) - p4.431) with 7 = 0 and g = oo of Proposition 14.51 respectively, 
using the nonlinear estimates (I4.44p - (|4.46l) . we deduce 

+ ||£>(t)||^oc + ||u(f)||^oo < (1+ (Ko+ 5A7(f)) . (4.69) 

We thus deduce from p4.68p - p4.69p that 

C-{t) < (1 + {Ko + 6M{t)) . (4.70) 

Now we substitute p4.70p into (I4.65P to obtain 

< e-^^Kl + [ + d^Af\t)) dr. 


(4.71) 


Finally, by the estimates p4.42p - (l4.43l) with 7 = 0 and g = 2 of Proposition 14.51 using the 
estimates p4(44p - p4.45l) with 7 = 0 and p4.67p . in view of (I4(64p . we obtain 

rt 


lll>(i)ll 


L 2 ^ 


(1 +1) Kq + 5 [ (1 + f —r) 2(’- 2 ) 2(l + r) ^ 'iJ\f{t)dT 

Jo 


< 

rsj 


(l + t)-^-5 {Ko + 


(4.72) 


and 

lk( 0 llL 2 < (1 + t) '^'^Ko + df (1 + 7 - r)“ 2 (F“ 2 )(i-g 7-)“‘’“4A7(t) dr 

Jo 

<(l + 7)-« (Ko + ~5M{t)). 

Note that we have used the fact (/ + | > 1 since p, r < 3/2 so that 

f (1 + 7 — t)~2(f“ 2)“2(1-)- r)“'’“4 dr < (1 + 7)“2(f“2)“2 

Jo 

and 

^ S/'llN 3 3/ll\ 

(1 + 7 — r)“2(F“2 -g t)~^“4 (;77 - < (1 -g t)~2\7~2). 


/ 


(4.73) 

(4.74) 

(4.75) 


By the definition p4.64p of A7(7), we deduce from (I4.70I) - (I4.73I) that 

Af{t)<Ko + ~6Af{t). (4.76) 

This implies 

A7(7) < Ko (4.77) 

since 6 is small by Proposition l2.ll This in turn together with the interpolation gives (ll.9l) - (ll.lll) 
for 6/5 ^ r < 3/2 by taking Cq = Kq. 

Case 2: 1 < r < 6/5. In this case, we define 

77(7) := ||p(7)||^2 + ||Vu( 7)||^2 + ||£'(7)||^oo + \\u{t)y^ , (4.78) 

+ (4.79) 
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and 


IC{t) := sup ({1 + t)^^ 2 (^U{t) + J{t))'). 

Osgr^t ^ ^ 


We take £ = 1 in (j4.59p to have, in view of ()4.78h - (j4.79p . 

<e-^^Kl+ f e-^^^-^^n\T)dT. 


(4.80) 


(4.81) 


We now estimate the time decay rates of 'H{t) by applying the linear decay estimates. By 
the estimates (|4.42|) with i = 0 and g = 2 of Proposition 14.51 and using the nonlinear estimates 
(I4.44p - (|4.45p . in view of (I4.78p - (|4.79p . we obtain 

I|£'WIIl2 < (1 +1) ^Ko+ f {l + t-T)~^(r- 2 )- 2 S['H + J){T)dT 

Jo 

dr. (4.82) 

<H\\q\\l^ + Mj,^)<6J. (4.83) 


Note that 


lL2 


V^-rU 


L 2 


+ \\u\ 


L 2 


V^-rU 


L 2 


Hence, we have 


I|£'(*)IIl 2 ^ (1 + ^) '^Ko+ [ (l+i-'r) 2 (-- 2 ) 2S{'H + J){t) dr 

Jo 

(1 + t) ^Kq + 6 f (1 + t — r)~2(““2)“2 (1 +-7-)“‘’“2/C(t) dr 

Jo 


< 

r\j 


< 


(l + t)-f-^ {Ko + 


(4.84) 


Here we have used the fact |(f“ 5 ) + ^>l since r < 6/5. Similarly, by the estimates p4.43p 
with i = \ and g = 2, and (I4.42l) - (l4.43p with £ = 0 and g = 00 of Proposition 14.51 respectively, 
using the nonlinear estimates p4.44l) - (|4.46p . we deduce 


l|Vn(t)||^2 + ||£»(t)||ioo + ||n(t)||^cx. < (1+ t) ^ 2 (^iTo + ^^(i)) • 
We thus deduce from (j4.84li - (l4.85h that 

Jiit) < (l + i)"^“^ {Ko + ~ 


Now we substitute (j4.86li into (I4.8ip to obtain 

rt 

e-^(*-")(l + T) 


(4.85) 

(4.86) 


f 


\K^ + 6‘^}C‘^it))dT. 


J\t) < e—+ 

By the definition (|4.80p of /C(t), we deduce from (|4.86p ~ (l4.87p that 

K:{t)<Ko + 51C{t). 

This implies, since 5 is small, 

m < K,. 

Finally, by the estimates (I4.43p with £ = 0 and g = 2 of Proposition 14.51 using the estimates 
(j4.44h - (l4.45p with £ = 0 and (I4.83p . in view of M.SOp . by (I4.89p . we obtain 

U[t )\\^2 + Ko + ~5 

0 

<Ko(l + t)-«. (4.90) 


(4.87) 

(4.88) 

(4.89) 


(1 +1) + S f (l + t — r) 2(7- 2) (^ir) ^ 4/C(t)dr 

Jo 
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Note that (|4.89p implies 


< ^o(l +t) ^ 2- 


(4.91) 


That is, we have proved (jl.9p for i = 0 and (ll.lOp for 0 ^ ^ ^ 1 by the interpolation. To 
prove the remaining decay estimates in (jl.9p - (ll.llh . we may now employ the arguments used 
in Case 1. Indeed, since 0<3 — 3/r<l and 4 — 3/r > 1, by the interpolation, we deduce from 
that 


I|£'('^)IIl 2 ru{T) + ||u(r )||^2 ?u(t) 


L2 


<K^{1 + t) 2 (l + r) ^ 2^ +K^{1 + t) ^(1 + r) ^ 2 




(4.92) 


Here we have used the fact 2 C + l/ 2 ^(( + 3/4 since p, r < 3/2. So by replacing the estimates 
(I4.67P by the estimates p4.92l) and then reproducing the arguments of Case 1, we may derive 

M{t)<Ko + Kl (4.93) 


where AA(t) is dehned by (|4.64l) . This in turn together with the interpolation gives (ll.9p - (ll.lip 
for 1 < r < 6/5 by taking Cq = Kq + Kq. 

Now in view of these two cases, the proof of Theorem 11.31 is completed. □ 


Appendix A. Analytic tools 

We recall the Sobolev interpolation of the Gagliardo-Nirenberg inequality. 

Lemma A.l. Let 2 ^ p ^ 00 and a, /3 ,7 G M. Then we have 

I|v“/|L. < v^/ (A.l) 

Here O^0^1 (if p = 00 , then we require that 0 < 0 < 1) and a satisfy 

a + 3 Q - = /3(1 - 6 ») + 70. 

Proof. For the case 2 ^ p < 00 , we refer to Lemma 2.4 in [5]; for the case p = 00 , we refer to 
Exercise 6.1.2 in [4]. □ 

We then recall the following commutator and product estimates: 

Lemma A.2. Let I ^ 0 and define the commutator 

V',p 


Then we have 




h 


< 


LPO 


llVpIl, 


In addition, we have that for I ^ 0, 
V\gh) 


LPO 


~ llffll LPi 


Here po,p 2 ,P 3 G (l,oo) and 


V\gh)-gyh. 


(A. 2 ) 

+ Vg 

1 ^ LP4 • 

(A.3) 

LP2 

LPO 

+ V^p 

^ LP4 • 

(A.4) 

LP2 LPO 

1 1 1 

P2 P3 PA ' 




Proof. We refer to Lemma 3.1 in [7]. 


□ 


Lastly, we record the estimates of the remainder TZ defined by (13.4p . 
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Lemma A.3. Let IZ he defined by (|3.4p . Then we have that for I 1, 


and 


V^TZ 


Vg 

Vtz 

^2 ~ 1 ,°° + 

Vg 


(A.5) 

(A. 6 ) 


Here “ < ” stands for “ ^ C” with the constant C depending on the function h, the upper and 
lower bounds of ps and || . 


Proof. We only prove (|A.5p . while (lA.Op can be proved similarly with minor modifications. We 
may view TZ as an operator over h, i.e., we define the operator TZ{f) of the smooth function /: 


Then TZ = TZ{h). 
Moreover, taking 


nf) ■■= [ 

jp. 

L 


f''{s){Q + Ps - s)ds 


= / f"{Q + Ps-T)TdT. 

Jo 

It is clear from the dehnition (IA.7h that 
7^(/) = 0{g^). 
the spatial derivative of ()A.7p yields 

V77(/) = f{ps)QVg + r f'"{g + ps - t)t dr (Vg + Vps) 

Jo 

^f"(ps)gVg + 7Z(f') (Vg + Vps). 


(A.7) 

(A. 8 ) 


(A.9) 


Hence, by Holder’s and Sobolev’s inequalities, we have 

I|V77(/)||^2 < \\nPs)QVQ\\L2 + \\nf) (V^? + Vp,)||^, 

< ||/"(p.)|Lco Ml- MqM + Ml- Ml^ (l|V^>||i3 + ||Vp,||^3) 

<Ml-\\^qM<^MqM- (A.IO) 


Since TZ = TZ{h), we deduce (IA.5h for / = 1. 

Now for / ^ 2, by the identity (IA.9h and the product estimates (IA.4h of Lemma IA.21 we 
obtain 

'VTZif) ^ = V^-^(V77(/)) 


< 

rs_/ 


L 2 

^'-HnPs)QVg) ^ + V-^TZ{f'){Vg + Vps)) 


L 2 


L 2 


L 2 


< V-\f{ps)) 

^3 + ||/''(Ps)||^oo 

V'-' {gVg) 

+ ll^(/')Loo 

V'l? 

,, + l|K(/')IL. 

Vps 

L3 

+ ( + W^Ps 

Ill») V'-'K(/') 

L 2 



< 


\ l 6 ||V£i||^oo + ll^'llioo V''g 


L 2 


+ 




1,6 


l|V^>||i3 


+ Ml- ^+Ml-Ml<^+ v^-W) 


L 2 


L 2 


< 5 


L 2 + 


V'^||^ 2 )+||^'''W) 


L 2 


By this recursive inequality pA.llh . we obtain that for I 2, 


(All) 


Vtz 


L 2 




e=i 


1,2 


1,2 


<'5(||Ve||^2 + ||v^f? 


L 2 


(A. 12 ) 
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Here in the last inequality we have used the inequality (lA.lOp for / = This proves (IA.5I) 

for I ^ 2, and the proof of the lemma is completed. □ 
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